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We test the possibility that the newly-discovered 126 GeV boson is a pseudoscalar by examining 
the correlations among the loop-induced pseudoscalar decay branching fractions to 77, ZZ* , Zf, 
and WW* final states in a model-independent way. These four decays are controlled by only two 
effective operators, so that the rates in Z'y and WW* are predicted now that the rates in 77 and 
ZZ*,Z^* — y 41 have been measured. We find that the pseudoscalar possibility is disfavored but 
not conclusively excluded. Experimental exclusion of the Z'y decay to well below a/asm ~ 170 
or conclusive observation of the WW* decay near the Standard Model rate would eliminate the 
pseudoscalar possibility. The Z7 exclusion should be possible using existing data. The only loophole 

£\| in our argument is the possibility that the 4£ signal comes from pseudoscalar decays to a pair of 

t-H new neutral gauge bosons with mass near the Z pole. 
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I. INTRODUCTION 

The ATLAS and CMS experiments at the CERN Large Hadron Collider (LHC) have discovered a new boson with 
mass of about 126 GeV [U |2] using about 5 fb _1 of integrated luminosity per experiment at each of 7 TeV and 
8 TeV. Decays into 77 and ZZ* — > At final states have been observed with significances of 4.1-4.5er [2] [3] and 3.2- 
3.4cr [UUj, respectively, at each experiment. These final states have good mass resolution of order 1 GeV. Decays into 
WW* — > tvtv have also been observed at the 1.6-2.8cr level [21 [5], albeit with much poorer mass resolution. Searches 
for decays into rr and bb final states have so far been inconclusive [H |H1 E] . The observed signal rates in all channels 
are consistent with the new boson being the Standard Model (SM) Higgs boson. 

A key question posed by this discovery is the determination of the spin and CP quantum numbers of the new boson. 
Observation of the two-photon final state excludes the possibility of spin 1 via the Landau- Yang theorem [8] Q Spin-2 
can be distinguished from spin-0 by examining the angular distributions of the final-state photons jTUl [TT] , the leptons 
in ZZ* — > At JTUl IT2~HT5] and WW* — > tvtv [TT] final states, and the angular correlations between the jets produced 
in vector boson fusion |16j . 

A spin-0 resonance can have CP = +1 (scalar, like the SM Higgs) or CP = —1 (pseudoscalar). The CP nature 
of a spin-0 resonance can be determined using the angular distributions of the leptons in ZZ* — > At [TH [T3J IT?rl2"2"] , 
the angular distributions of the jets produced in vector boson fusion [53] or gg — ► Xjj [51] . hadronic event shape 
observables [55], or the spin correlations in decays to rr [26J. All of these measurements require more integrated 
luminosity than has been used up to now for Higgs analyses. 

Several pseudoscalar interpretations of the 126 GeV excess, particularly in the 77 final state, have been put forward 
in the context of specific models [37H33]. A scalar-pseudoscalar admixture has also been considered in Ref. [54] . 

In this paper we attempt to constrain the pseudoscalar possibility in a model-independent way by examining the 
correlations among the 77, ZZ*, WW*, and Z7 decay modes. For a pseudoscalar, all of these decays arise from 
loop-induced effective couplings, which can be parameterized in the electroweak basis in terms of just two operator 
coefficients. A similar approach has been taken to constrain an electroweak-singlet CP-even scalar that has only 
loop-induced couplings to WW and ZZ in Ref. [35] . We improve upon the analysis in Ref. 35j by taking into account 
the very large contribution to the At final state coming from — > Zj* — > At. 

After fitting the ratio of the 77 and Al final states to the SM Higgs prediction, we find that the branching fraction 
for the pseudoscalar into WW* is dramatically suppressed compared to the SM Higgs and the branching fraction 
into Zj is dramatically enhanced. The suppression of the WW* mode disfavors the pseudoscalar at somewhat less 
than the 3<r level, based on the ATLAS excess in this channel [5J. The enhancement of the Zj mode is not large 
enough to allow the pseudoscalar to be excluded based on the LHC measurements of the continuum Zj production 
cross section; however, it is large enough that a resonance search should easily have the sensitivity to exclude the 
pseudoscalar possibility. 

We also consider possible loopholes in the argument for the large and easily-excludable enhancement of the Zj 
branching fraction. We conclude that the only way out is if the observed At final state comes from decays of the 
pseudoscalar into a pair of new neutral gauge bosons with mass near the Z pole, which then decay to t + £~ . The 
strong suppression of the WW* branching fraction survives as a distinguishing feature in this case. 

This paper is organized as follows. In Sec. [TTJwe define the effective operator parameterization for the pseudoscalar 
couplings to gauge bosons. In Sec. |III| we determine the parameter value needed to match the observed rates of the 
126 GeV boson into 77 and At final states, taking care to properly include contributions from cf> — > Zj* — > At. In 
Sec.[lV]we compute the resulting predictions for <j) — > WW* and <p — > Zj and discuss how they can be used to confirm 
or exclude the pseudoscalar possibility. In Sec. [V] we consider potential loopholes in our analysis. We conclude in 
Sec. ETJ 

II. PSEUDOSCALAR COUPLINGS 

A pseudoscalar has no gauge-invariant renormalizable couplings to SM gauge fields. Instead, we model the interac- 
tion of the pseudoscalar with the SM gauge fields using the following effective Lagrangian: 

B^B^ + bW^W^] . (1) 

Here v = 246 GeV is the SM Higgs vacuum expectation value. The dual field strength tensors are defined by 
qo-hv — ^ u P a G a p<y , and similarly for the other field strength tensors. The above Lagrangian allows us to parameterize 



1 A spin-1 boson could still be accommodated in the scenario that the 77 signal comes from decays to two very light intermediary particles 
which each decay to a pair of tightly collimated photons [9] . 



3 





0.0001 r 




1 e-ut> r 




1 e-uo r 




1e-07 r 


> 




(Ge 


1e-08 r 








1e-09 r 




1e-10 r 




1e-11 - 




1e-12 - 







yy\ 




^^^^ 








1 1 1 1 1 1 







-10 




b 



10 



FIG. 1. Theoretical prediction for the partial widths of cj> to 77 (blue), Z7 (red), WW* (green), and ZZ* (magenta) as a 
function of b, for M$ — 126 GeV and a = 1. The normalization on the y-axis can be changed by varying a; all the partial 
widths shown scale with a 2 . The vertical dotted line shows the value of b that we will use to match the SM rates in 77 and 4£. 

the decays of tj> to 77, W + W~, ZZ, Z-y, and gg. The amplitudes for these processes can be easily derived from 
Eq. At leading order, their dependence on the parameter 6 is: 

M (<f> -> WW) oc 6 
M{(f>^ ZZ) oc {bca^Qw + sin 2 9 w ) 
M(4> -> 77) ex (cos 2 ^ + fesin 2 ^) 

Z7) cx (6- 1), (2) 

where 0^ is the weak mixing angle. 

It is readily apparent from the above equation that the decay amplitude to a particular gauge boson pair can be 
adjusted by tuning the parameter b. In particular, each amplitude vanishes at a different value of b: 

M((j}^ WW) = when 6 = 

M {(j) ZZ) = when b = - tan 2 9 W 

M((j> -> 77) = when b = — cot 2 9 W 

M (0 -> Z7) = when 6=1. 



(3) 



We illustrate this behavior by plotting the partial widths of </> to Z7, 77, ZZ* , and M^W 7 * as a function of 6 in Fig. [TJ 
for a 4> mass of 126 GeV. 

Notice in particular that the partial widths to WW* and ZZ* are in general much smaller than those to 77 or Z7, 
except for 6 values near the zeros of the 77 or Zj amplitudes. This is due to the kinematic suppression of these decays 
for (j> masses below the WW and ZZ thresholds. This kinematic suppression will force us to choose 6 very close to the 
zero of the 77 amplitude in order to match the rates in 77 and 4£ experimentally observed for the 126 GeV boson. 



III. MATCHING THE OBSERVED RATES 



We now determine the value of 6 required to fit the observed rates in 77 and At final states of the 126 GeV boson. 
This can be done using only the ratio of rates (the overall production rate for gg — > can be adjusted by varying c). 
The fit is complicated by the fact that 4> Zj* 4^ contributes substantially to the At final state selected by the 
LHC experiments. 
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M 12 (GeV) M 34 (GeV) 

FIG. 2. Invariant mass distributions M12 (left) and M34 (right) for the 4e final state event selection (see text for definitions), 
after application of the A_R cut. We show separately the contributions from <f> — > Z"/* — > 4c (red) and (f> — > ZZ* — > 4e (blue, 
multiplied by 1000). The vertical lines indicate the lower cuts on M12 at 50 GeV and on M34 at 20.5 GeV. 



We proceed as follows. In the SM, the ratio of Higgs partial widths into ZZ* versus 77 for a Higgs mass of 126 GeV 
is predicted to be [36 PI 



-,s.\j _ -> ZZ *) 



This can be rewritten in terms of a partial width after cuts in the 4e final state as 

SM T(H -> ZZ* -» 4e) 1 y r c (ff -> ZZ* -» 4e) 

"ss/rr r -(JT -+ ZZ* 4e) [ B R(Z ee)] 2 r(if ^ 77) 



where r c refers to the width after cuts (defined in terms of the total signal rate after cuts; the boson production cross 
section and total width factors cancel in the ratio with T(H — > 77)). The first term in Eq. ([5| is determined by a 
Monte Carlo simulation using CalcHEP [37] , the second is computed at leading order to match the signal calculation, 
and the third is what has been measured by ATLAS and CMS. 

We match the observed rates by requiring that the ratio of rates to A£ and 77 for the pseudoscalar be the same as 
that of the SM Higgs; i.e., in Eq. jjj) we make the replacement, 

T C (H ZZ* -> 4e) T c (0 -> 4e) 

t(h 77) ^ r(0 ^ 77) ■ ( ^ 

To determine the acceptance for H —> ZZ* —> 4e and the partial width after cuts for 4> — > 4e, we apply a subset of 
the ATLAS selection cuts g]: 

• Each pair of leptons must be separated by at least Ai? = 0.1. 

• The electron-positron pair with invariant mass Mi 2 closest to the Z mass must satisfy 50 GeV < M12 < 106 GeV. 

• The remaining electron-positron pair must have an invariant mass M 34 between 20.5 GeV and 115 GeV. This 
cut affects how much the process Z7* —> 4e contributes to the signal after cuts, as can be seen in Fig. [2] Note 
that the Z7* contribution after cuts is about 2400 times the size of the true ZZ* contribution, in extreme 
contrast to the SM Higgs. 

Putting everything together and adjusting the value of b until Eq. ([5| is satisfied, we obtain two solutions for b, 
corresponding to the two sides of the zero in the 77 rate shown in Fig. [T] For concreteness we take the side with 
smaller \b\, 

b= - cot 2 9 w (l + e), where e = -0.092. (7) 



2 The contribution of SM H — > Z7* to the Al final state is experimentally negligible because F(H — > Zy)/r(H — ► ZZ*) = 0.056; taking the 
photon off shell suppresses the Z-y* contribution even further. 
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Note that b must be tuned to within ten percent of — cot 2 8w- Because of the rapid variation of T(<fi —> 77) in this 
range of b, factor-of-two changes in the ratio T c (<f> — > 4e)/T(4> — > 77) result in only small absolute changes in the value 
of b. Our results in the next section for the rates for — » WW* and <f> ~ > Z7 do not depend sensitively on such small 
variations in b. This justifies the several assumptions that we have made in the above analysis: 

• We fit the value of b to the SM prediction for T(H — > ZZ*)/T(H — > 77), rather than to the measured ratio. The 
measured ratio is smaller than the SM prediction, but consistent within the experimental uncertainties (ATLAS 
and CMS both measure a rate in the 77 channel somewhat larger than the SM expectation) . Fitting to the SM 
prediction rather than the measurements yields a slightly different value of b but will not change our overall 
conclusions. 

• We assume that the experimental acceptance for the pseudoscalar <fi — > 77 signal is the same as that for the 
SM H — > 77 signal. This allows us to write T(H — > 77) in the denominator of Eq. Q and T(cf> — > 77) in 
the denominator of Eq. ^ without including selection efficiencies. Because the decay kinematics are identical, 
the difference in efficiencies in this case can come only from the difference in the pt distribution of the initially 
produced boson and the kinematics of the accompanying jets. A proper evaluation of the difference in efficiencies 
would require a full simulation of Higgs and pseudoscalar production, which is beyond the scope of our analysis. 
In any case, this effect is small enough that it will not change our conclusions. 

• We do not apply pr or rapidity cuts to the leptons in the 4e final states. This will affect our results only to the 
extent that H and <j) are produced with different transverse momentum distributions. Again, this effect is small 
enough that it will not change our conclusions. 



IV. PREDICTIONS FOR WW* AND Z 1 FINAL STATES 



With b in hand, we can read off the predictions for the pseudoscalar decay widths to WW* and Z7. We find, 

< w/ ^->ww* ) = 
zi/11 r(0^7 7 ) 1 ' 

For comparison, the corresponding ratios for the 126 GeV SM Higgs are [33] 

R S ww hl - 101, R S ^ hl * 0.711. (9) 

The smallness of the rate for <fi —> WW* , 440 times smaller than the SM H — > WW* rate, disfavors the pseudoscalar 
insofar as ATLAS and CMS have observed weak signals in this channel (1.6<r from CMS and 2.8a from ATLAS [H[S]). 
In particular, the ATLAS analysis of 8 TeV data, which provides the largest contribution to their overall significance, 
so far includes only the e/i final state. This is interesting because Z7* — > eevv, \i[wv will spill in to the same-flavor 
event selections, thereby boosting the rates in these channels, but it cannot contaminate the ev\iv sample. 

The most direct probe, however, is the very large enhancement of the rate for pp — > <j) —¥ Z7, 170 times larger than 
the SM pp — >• H —> Zj rate. This translates [3S] into a cross section for pp —> <fi —> Zj of 4.8 pb (6.0 pb) at 7 TeV 
(8 TeV) . This cross section is too small to be constrained by direct comparison with the continuum Z-f cross section, 
measured by ATLAS and CMS to be about 290 pb (up to 10% uncertainty) in the fiducial region of Ej, > 10 GeV, 
AR e ^ > 0.7, M u > 50 GeV [3JHD]. 

On the other hand, <j) — > Z^ — > M'y produces an invariant mass peak. A phenomenological study of this channel has 
been done for the SM Higgs in Ref. (41], which found a 95% CL exclusion reach of ct/ctsm — 4 for a 126 GeV Higgs 
with 20 fb _1 at 8 TeV after an optimized cut on a multivariate discriminant (including only statistical uncertainties). 
Scaling to the luminosity of ~ 5 fb _1 already used for Higgs analyses at 8 TeV, this should translate into an idealized 
exclusion reach of about ct/ctsm — 8- The pseudoscalar prediction of ct/ctsm — 170 in the Zj channel should therefore 
be easy to either discover or exclude in current data, even with a non-optimized analysis. 

Finally we comment on the total width of the pseudoscalar. So long as BR(0 — > gg) -C 1, the total width is 
dominated by Z7 and we find BR(</> — > 77) ~ 1/121 = 0.83% ~ 3.6BR(iJ — > 77). In this case the signal rate in 
pp — » 4> — > 77 matches the SM Higgs prediction for a(gg — > <f>) ~ 0.3 cr{gg — > H). This can be used to fix c. The overall 
scaling factor a on the effective electroweak gauge couplings can then be adjusted to maintain BR(</> — > gg) -C 10 



The parameters c and a can of course be increased together to maintain a constant value of cr(pp — ¥ (ft) X BR(</> — > 77) while increasing 
the total width of cj>. 
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Insofar as BK(H — > gg) = 0.085 <C 1 in the SM [3B], the total width of the pseudoscalar can therefore be comparable 
to the SM Higgs total width of about 4.2 MeV [36], i.e., well below the detector resolution. 

V. POTENTIAL LOOPHOLES 

We can think of only two potential loopholes in our prediction for the easily-excludable large enhancement of the 
Zj rate for a pseudoscalar. Both involve the breakdown of the effective Lagrangian description of Eq. ([!]). 

A. Dependence of the effective couplings on gauge boson masses 

Starting from a renormalizable theory, the effective couplings of the pseudoscalar to a pair of gauge bosons in 
Eq. (JlJ are generated by loops of fermions that couple to the pseudoscalar and carry the appropriate gauge charge. In 
the limit that the fermions are much heavier than the pseudoscalar, W, and Z, the effective Lagrangian is recovered. 
However, if the fermions are light enough, the loop amplitude will depend on the invariant masses of the external 
gauge bosons. In particular, the coefficients a of the hypercharge effective coupling and ab of the SU(2) i effective 
coupling could take different values for the Zj final state than for the 77 final state. This has the effect of enhancing 
the contribution of any particular fermion to Z7 compared to its contribution to 77, as follows. 

The fermion loop function for on-shcll cf> — > 77 is given by |42j 



and for on-shcll <p — > Zj by 



Ff /2 {r) = -2t/(t) -2 for m/ ->■ 00, (10) 



7 2 A ) = - o7^~~~T\ lf(r) - /(A)] \ for m f -> 00, (11) 



2(r — A) J WJ 2 



where r = 4ml /M?, A = 4m*/M§, and 



f(x) = i L sm wv-;j > (12 ) 



'(V^)] 2 r>l, 



, ~I [MW 7 ?-) - T < 1, 

with r)± = 1 ± yl — x. Forming the ratio (normalized to 1 at m/ — > 00 



, h(r, A) f 1.087 for m f = 100 GeV 



4 .F/ 2 (t) ( 1.018 for m / = 200 GeV ^ 13 ^ 

This effect can be used to shift the effective b value in the Z-f amplitude upward (for a light fermion carrying only 
weak isospin) or downward (for a light fermion carrying only hypercharge). We see that the shift in this effective b 
value is constrained to be less than 9% in either direction (we take 100 GeV as a lower bound on the mass of such a 
new fermion due to its non-observation at LEP-2). Such a shift in the effective value of b is too small to significantly 
change the prediction for T((f> — > Zj), so it does not provide a loophole in the exclusion. 

B. Decays of <f> to weak-scale Z' boson pairs 

Our prediction of the huge enhancement in the Z"f rate for the pseudoscalar rests on the requirement to tune the 
077 effective coupling to be close to zero in order to match the observed ratio of rates in the 77 and 4£ final states. But 
what if the observed 4£ final state does not come from the effective vertices in Eq. ([I])? In this case, b is unconstrained 
and we can have T(<f> — > Zj) ~ T{<p —> 77) as it is for the SM Higgs. 

This can be achieved if the observed events in the 4£ final state come from cf> — > Z'Z' — > 4£, for a new neutral gauge 
boson Z' with mass near the Z pole (in order to approximate the kinematics of the SM H — > ZZ* — > 4£ decay) . Such 
a Z' under the Z pole is in fact slightly favored by LEP measurements [13] . 

In this case, a conclusive detection of the WW* final state at a rate consistent with the SM Higgs would still 
exclude the pseudoscalar possibility, which predicts T(<fi — > WW*) to be more than two orders of magnitude smaller 
than T((j) — ¥ 77) for generic b values not too close to the zero in the 77 rate. 
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VI. CONCLUSIONS 

In this paper we examine the possibility that the newly-discovered 126 GeV boson is a pseudoscalar by parame- 
terizing the couplings to gauge bosons in terms of an effective Lagrangian. In this framework, the four decay modes 
to 77, ZZ* , Z7, and WW* are controlled by only two effective couplings, i.e., the coefficients of the (pB^B^ 

and 4>W^W lfJ,l/ effective operators. The ratio of these coefficients can be fixed in terms of the ratio of the highest- 
significance experimentally-observed rates into 77 and A£, leading to predictions for the rates in WW* and Zj. 

An important part of our analysis is the realistic treatment of the Zj* — > 4£ contribution to the 4£ final state 
after experimental cuts. In fact, we find that Zj* dominates over ZZ* by a factor of ~ 2400 after cuts. Taking this 
into account, the pseudoscalar model predicts a rate ct/ctsm — 170 in the Zj final state relative to the SM Higgs. 
While this enhancement is not large enough to show up against the SM Z7 total rate, it should be easily visible (or 
excludable) through a Z-f resonance search using only the existing LHC data. This prediction can be evaded if the 
observation in the At final state is due to exotic decays of the pseudoscalar to weak-scale Z' pairs. 

We also find that the decay rate to WW* for the pseudoscalar is dramatically suppressed compared to that for the 
SM Higgs, yielding <t/<t sm ~ 1/440 for this channel in the ev\xv final state, which is not contaminated by Z-f* decays. 
Conclusive observation of the WW* final state with a rate near the SM Higgs prediction would thus eliminate the 
pseudoscalar possibility. This, however, will require additional integrated luminosity beyond that used already in the 
Higgs analyses. The prediction of the suppression of WW* is robust against new contributions to the 4£ final state 
from exotic decays of the pseudoscalar to weak-scale Z' pairs. 
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